x=[0:2:10]';
y=[0 2.9 14.8 39.6 74.3 119]';

A=[ones(size(x)) x x.^2 x.^3 x.^4 x.^5];
a=A\y

A=[ones(size(x)) x x.^2 x.^3 ];
a=A\y
A*a-y

% x=[1 2 3];
% y=[12 15 16];
% plot(x,y,'o');
% hold on
% t=0:0.01:4;
% plot(t,a(1)+a(2)*t+a(3)*t.^2,'r')
% xlabel('t')
% ylabel('p(t)')
% axis([0 4 7 17])